
National Aeronautics and Space Administration 


CAUCHY DRAG ESTIMATION IjlMK 
EARTH ORBITERS 

* * ' . 

• * J. Russell Carpenter 


I ^ Alinda K. Mashiku 


NAVIGATION & MISSION DESIGN BRANCH 

NASA GSFC 


www.nasa.gov 



GPS-based Ow 
Sometimes Inf 


Initial JSpOC 
prediction 
~2x closer to 
final 

consensus at 
Time of 
closest 
approach 
(TCA) 


Initial 
predictions 
similar; GPS 
moving 
toward 
consensus 
sooner 


Case 1 



0 

185 186 187 188 189 190 191 192 193 194 195 

Day of Year - CARA Report Date 


Case 3 



Day of Year - CARA Report Date 



Initial JSpOC 
prediction 
~3x closer to 
final 

consensus at 
TCA 


Initial JSpOC 
prediction 
within 500 of 
final 

consensus at 
TCA 


Ref: M.A. Vavrina, C.P. Newman, S.E. Slojkowski and J.R. Carpenter, “Improving Fermi Orbit Determination and Prediction in an Uncertain Atmospheric Drag Environment” 
Proceedings of the 24 th International Symposium on Space Flight Dynamics, www.issfd.org , 2014 



NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


2 




Max 3-Day Position Prediction Error Magnitude (m) 


Tuning GPS E 

m 

arginal Impr 

ovement 

in Prediction 1 


\ to Density V 

ariation 


3000 


2500 


2000 


1500 


1000 


500 



0 *** ^ V r . T 

v to v v q> /\> > x a? C+ n> ^ o> -> ^ A x ^ A> of- _%> *f A x _of jV <S> ^ ^ ^ *- 

Solution Epoch 

Newman, S.E. Slojkowski and J.R. Carpenter, “Improving Fermi Orbit Determination and Prediction in an Uncertain Atmospheric Drag Environment” 
International Symposium on Space Flight Dynamics, www.issfd.org , 2014 



Ref: M.A. Vavrina, C.P. Newman, S.E. Slojkowsk 
Proceedings of the 24 th International Symposium 


? >H? ^ ^ ^ ^ ^ ^ ^ 

<f <f ^ <f # # # <f # ^ 

'oj'ofof't^opsvtvcvc^ o^ 


1000 


950 


900 


850 


800 


750 


NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


3 


Exospheric Temperature (K) 


Solar Flux Variati 
Activity Intervals 


Data Set: Observed F 10 7 Flux when ISN >75 


Variations in Data Set: Time differences of 

/ n\ : i . . /\l i n a _i _ . . _ . . _*\ . r\±. i n a _i _ . . _ . . 


F 1Q7 forlSN >75 



50 - 

Q | 1 1 J . i—l 1 J 1 1 l—i 

1955 1960 1965 1970 1975 1980 1985 1990 1995 2000 2005 2010 2015 



Time differences of Double Differences 




NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


4 


Drag Residua 



aussian 


0.002 


c<r 0.0015 

1 0.001 
CO 

o 0.0005 
ro 

0 0 

0 

o 

ra-0.0005 

(0 

3 

2 - 0.001 

o 

-0.0015 


- 0.002 

11-5-04 11-6-04 11-7-04 11-8-04 11-9-0411-10-0411-11-04 




NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


5 


Probability Density 



Fit PHFq f nr F Rata 




Stable Distribution Fits 


Data 

Concen- 

tration* 

Asym- 

metry 

Scale 

^10.7 

1.38 

0.02 

0.13 

a p 

1.03 

0.0014 

3.46 

CHAMP 

drag 

1.34 

0.6121 

0.0001 

*Gaussian 

= 2.0, Cauchy 

= 1.0 



NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


6 



The Idan-Speyer 


mator (ISCE) 


■ Given a linear scalar system 


%k + 1 = 4>k^k + w k 

Vk = H k x k + v k 


■ Posterior PDF 


k+2 


Px*|Y k (a*|Y*) = 




^(0ife|fc) 2 + w (*)fc|fc 


■ With Cauchy inputs and initial 
condition 


Pxn(^0) 


a/ 7r 

(xo - xo ) 2 + a 2 

(3/ir 


Pw '■(“'*) w 2 + /? 2 

Pv t ( ''/;) = 2~T~~ 2 

+ 7 


Has Finite Moments 


-oo 


£fc|fc = E[xfc|Y fc ] = / &p Xfc |Y fc (£fcl Y fc) d f 


-OO 


oWfelfe (o’Wfclfc) 2 - + &Wfc|ife< 7 Wfc|fc 

= 7r !L 

i—\ 


CjJ 


^})k\k 



NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


7 



11 


Truth 

Estimate 

— — - Est.± 1 o 
PDF(± 6a) 



vJ\e& s 


8 




The ISCE May 
Density Estim; 



If I % i 


ie EKF for 




NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


9 



Consider States: c k \ k = x k \ k = ISCE mean 

Pss klt = (l~K s [H s H e \)P, 


S&k\k — 1 sP- cPcS k | k — 1 


Psc k , k — P CSklk _ 1 — {i Ks [h s h c ]) p SCk | fc _! iir s iy c p CCfe|A ,_ 1 


Pcc klk =Px k | fc = ISCE variance 


P P 

ss fc + l|A: sc /c+l|fe 

P P 

C5 fc + l|/c CC /c + l|/c 


'fclfc 

Sfc+llfc 

.£fc+llfc. 


= ^te+lp/c) 


r f i 

tk 

K L a ( r ). 

p 

ss k\k 

p 

sc k\k 

P 

cs k\k 

P 

cc k\k 


dr 





NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


10 


Error in Actual In-track Radial 

Density Density Perturb. OD Error OD Error 


Definitive OD is Si 
Estimation is Sup 


Dashed lines = ±3o formal error 


Baseline EKF 



J l l l l L 

t 1 1 1 1 r 



j i i i i L 



EKF Disciplined by ISCE 



t 1 1 1 1 1 r 



J i i i i i L 


1 1 1 I - 

1 1 1 

1 1 1 1 

1 1 1 


1 1 1 1 

i i 

T 

i , 

1 K 1 ■ 1 r * M 

i i i i_ 

1 v‘ ln * T *’ ^ ^ T ’ { 

i i 

j 



NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


11 


llft'fT 




Error in Actual In-track Radial 

Density Density Perturb. OD Error OD Error 


Predictive OD: 
Robust to Uno 




and 


sions 


Baseline EKF 


Dashed lines = ±3o formal error 
* Observed during definitive span 

EKF Disciplined by ISCE 




NAVIGATION & MISSION DESIGN BRANCH, CODE 595 

NASA GSFC 


12 


Conclusions 


■Space weather data show heavy-tailed characteristics 
that are better modeled by Cauchy than Gaussian 

■Cauchy estimator (ISCE) may be embedded in EKF, 
using Schmidt-Kalman consider framework, for 
density estimation 

■Definitive OD performance indistinguishable from EKF 

■Predictive OD performance superior to EKF 
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